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Abstract 

We have searched for the r lepton flavor violating decays r~ — > i~Kg {I = e or /i), using a data 
sample of 281 fb _1 collected with the Belle detector at the KEKB e + e~ asymmetric-energy collider. 
No evidence for a signal was found in either of the decay modes, and we set the following upper 
limits for the branching fractions: B(t~ —> e~K$) < 5.6 x 10 -8 and B(t~ — > /U~Kg) < 4.9 x 1CT 8 
at the 90% confidence level. These results are improvements by factors of 16 and 19, respectively, 
compared with previously published limits from CLEO. 

PACS numbers: 
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INTRODUCTION 



Lepton flavor violation (LFV) is allowed in many extensions of the Standard Model (SM), 
such as Supersymmetry (SUSY) and leptoquark models. Lepton flavor violating decays with 
Kg mesons occur in models with either heavy singlet Dirac neutrinos pj, R— parity violation 
in SUSY or dimension-six effective fermionic operators that induce r — \i mixing 
Experiments at the I?-factories allow searches for such decays with a very high sensitivity. 
The best upper limits of B{r~ -> e~K%) < 9.1 x 1(T 7 and B{t~ -> < 9.5 x 10" 7 at 

the 90% confidence level were obtained in the CLEO experiment using 13.9 fb" 1 of data 

In this paper, we report on a search for LFV r decays r~ — > £~K§ {£ = e or /x)[f] with 
a data sample of 281 fb _1 collected at the T(4S f ) resonance and 60 MeV below it with the 
Belle detector at the KEKB e + e~ asymmetric-energy collider 0. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 

Particle identification is very important in this measurement. Particle identification likeli- 
hood variables are based on the ratio of the energy deposited in the ECL to the momentum 
measured in the SVD and CDC, the shower shape in the ECL, the particle range in the 
KLM, the hit information from the ACC, the measured dE/dX in the CDC and the parti- 
cle's time-of-flight from the TOF. For lepton identification, we use a likelihood ratio based on 
the electron probability V(e) 0] and muon probability V(fJi) (H determined by the responses 
of the appropriate subdetectors. 

For Monte Carlo (MC) studies, the following programs have been used to generate back- 
ground events: KORALB/TAUOLA B for r + r~, QQ [u} for BB and continuum, BH- 
LUMI jll| for Bhabha events, KKMC Q for e+e" -> /i+fi' and AAFH for two-photon 
processes. Signal MC is generated by KORALB/TAUOLA. Signal r decays are two-body 
and assumed to have a uniform angular distribution in the r lepton's rest frame. The Belle 
detector response is simulated by a GEANT 3 Q based program. All kinematic vari- 
ables are calculated in the laboratory frame unless stated otherwise. In particular, variables 
calculated in the e + e~ center-of-mass (CM) frame are indicated by the superscript "CM" . 



DATA ANALYSIS 

We search for r + r~ events in which one r (signal side) decays into £K® where I is e or 
/i and K$ decays into 7r + 7r~, while the other r (tag side) decays into one charged particle 
of opposite sign to the lepton with any number of additional photons and neutrinos. Thus, 
the experimental signature is: 

{t~ -> £-{= e~ or fi~) + ifgH ^ + ^)} + {r + -> (a track)+ + (n^ AG > 0) + X(missing)}. 



[|] Unless otherwise stated, charge conjugate decays are implied throughout this paper. 
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All charged tracks and photons are required to be reconstructed within the fiducial volume 
defined by —0.866 < cos# < 0.956, where 9 is the polar angle with respect to the direction 
opposite to the e + beam. Charged tracks should have momentum transverse to the e + beam 
p t > 0.1 GeV/c and photons should have energies E 1 > 0.1 GeV. 

Signal events are isolated with following selections. We first demand that the four tracks 
have zero net charge. The magnitude of the thrust jT^j is required to be larger than 0.9 
to suppress the qq continuum background. The event must have a 3-1 prong configuration 
relative to the plane perpendicular to the thrust axis. The K$ is reconstructed from two 
oppositely-charged pions that have an invariant mass within 0.482 GeV/c 2 < M n + W - < 0.514 
GeV/c 2 . The ir + iT~ vertex is required to be displaced from the interaction point (IP) in the 
direction of the pion pair momentum . In order to avoid fake K$ candidates from photon 
conversions (i.e. 7 — > e + e~), the invariant mass reconstructed by assigning the electron mass 
to the tracks, is required to be greater than 0.2 GeV/c 2 . We apply the lepton identification 
requirements to each track except for the two tracks that are part of the K§ candidate on 
the signal side. The electron and muon identification criteria are Vie) > 0.9 with p > 0.3 
GeV/c and V(fi) > 0.9 with p > 0.6 GeV/c, respectively. After the event selection described 
above, most of the remaining background comes from t + t~ and continuum events including 
a single K$ meson. 

To ensure that the missing particles are neutrinos rather than photons or charged particles 
that are outside the detector acceptance, we impose additional requirements on the missing 
momentum vector, p m { SS , calculated by subtracting the vector sum of the momenta of all 
tracks and photons from the sum of the e + and e~ beam momenta. We require that the 
magnitude of p miss be greater than 0.4 GeV/c and that it point into the fiducial volume 
of the detector: —0.866 < cos^^ < 0.956, as shown for the r~ — > fA~K$ mode in Fig. [T] 
(a) and (b). The total visible energy in the CM frame E™, is defined as the sum of the 
energies of the K$ candidate, the lepton, the tag-side track (with pion mass assumed) and 
all photon candidates. We require it to satisfy 5.29 GeV < E™ < 10.0 GeV (see Fig. [T] 

(c) ). Since neutrinos are emitted only on the tag-side, the direction of p m i SS should lie within 
the tag-side of the event. The cosine of the opening angle between p miss and the tag-side 
track in the CM system, cos ^tag— miss * s therefore required to be greater than (see Fig. Q 

(d) ). To suppress the continuum background, the following requirements on the number of 
the photon candidates on the signal and tag side are imposed: n^ IG < 1 and n^ AG < 2, 
respectively. 

Finally, the correlation between the reconstructed momentum of the £K^ system, peK s , 
and the cosine of the opening angle between the lepton and K<$, cos 6ik s is employed to 
exclude background from generic r + r~ and continuum: cos6(.k s < 0.14xlog(p£A' s — 2.7)+0.7, 
where Pik s is in GeV/c (see Fig. |2J). While this condition retains 99% of the signal events, 
99% of the generic r + r~ and 84% of the uds continuum background are removed. The signal 
detection efficiencies for the r~ —> e~K$ and r~ — > fi~K^ modes are 15.0% and 16.2% after 
all criteria applied, respectively. 

RESULTS 

Signal candidates are examined in the two-dimensional space of the £~K§ invariant mass, 
M inv , and the difference of their energy from the beam energy in the CM system, AE. A 
signal event should have M inv close to the r-lepton mass and AE close to 0. For both modes, 
the M inv and AE resolutions are parameterized from the MC distributions around the peak 
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FIG. 1: Kinematical distributions used in the event selection after K§ mass and muon identification 
requirements: (a) the momentum of the missing particle; (b) the polar angle of the missing particle; 
(c) the total visible energy in the CM frame; (d) the opening angle between the missing particle and 
tag-side track in the CM frame. The signal MC distributions are indicated by the filled histograms, 
the total background including r + r~ and qq is shown by the open histogram, and closed circles 
are data. While the signal MC (r~ — > fi'K^) distribution is normalized arbitrarily, the data and 
background MC are normalized to the same luminosity. Selected regions are indicated by arrows 
from the marked cut boundaries. 



with bifurcated Gaussian shapes to account for initial state radiation. These have widths 
ajgf™ = 6.2/7.4 MeV/c 2 and a^| h/low = 20/26 MeV, for the r~ -> e~K% mode and, 
^wgh/iow = 6 ^5 g Me v/ C 2 and a^p 1 ™ = 19/23 MeV, respectively for the r~ -> ^Kg 
mode, where the "high/low" superscript indicates the higher/lower side of the peak. 

We blind a region of ±5o"M inv around the r mass in M inv and a region of —0.5 GeV < 
AE < 0.5 GeV so as not to bias our choice of selection criteria. Figure |3] shows scatter-plots 
for data and signal MC samples distributed over ±15er in the M inv — AE plane. Most of 
the surviving background events in both modes come from — > tt ± K^ and D —>■ i^vK^ 
decays. The remaining continuum backgrounds in the r~ — > [i~K® mode are combinations 
of a true K$ meson and a fake lepton. 

To optimize the sensitivity for the search, we use an elliptically shaped signal region, which 
has major and minor axes that correspond to ±5cr in the MC resolution for the M inv — AE 
plane. This region is shown in Fig. |3] (a) and (b). Signal efficiencies for this region after all 
requirements are 11.8% for t~ — > e~K§ and 13.5% for r _ — > fi~K§, respectively. 

As there are few remaining MC background events in the signal ellipse, we estimate the 
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FIG. 2: Scatter-plots of (a) signal MC (r" -> n'Kg), (b) continuum MC, (c) generic r+r" MC 
events and (d) data on the PtK s vs cos 0£k s plane. Selected regions are indicated by curves with 
arrows. 



background contribution using the M inv sideband regions. This is achieved by extrapolating 
to the signal region under the assumption that the background distribution is flat in M inv . 
We find the expected background in the ellipse to be 0.2 ± 0.2 events for each of the two 
modes. We open the blinded region and find no data events for the t~ — > e~K§ and 
t~ — > f^~Kg modes (see Fig. |3] (a) and (b)). Since no statistically significant excess of 
data over the expected background in the signal region is observed, we apply a frequentist 
approach to calculate upper limits for the signal yields [T3|. The resulting limits for the 
signal yields at 90% confidence level, S90, are 2.23 events in both modes. The upper limits 
on the branching fraction before the inclusion of systematic uncertainties are then calculated 
as 

(1) 



B(r 



2eN tt B{K$ -> 71+71-) 



where N TT 



250 x 10 6 and B(Kg 
) < 5.5 x 10~ 8 and B{r 



71 1 71 



= 0.6895 1 
< 4.8 x 10 



The resulting values are 



The dominant systematic uncertainties on the detection sensitivity: 2eN TT B{K. 



7X + 7i~ 



come from Kg reconstruction and tracking efficiencies. These are 4.5% and 4.0%, 
respectively, for both modes. Other sources of the systematic uncertainties are: the trigger 
efficiency (0.5%), lepton identification (2.0%), MC statistics (0.3%), and luminosity (1.4%). 
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FIG. 3: Scatter-plots of data in the Mi nv - AE plane: (a) and (b) correspond to the ±15u area 
for the r~ — > e~Ko i and r~ — > ^"-Kg modes, respectively. The elliptical signal region shown by a 
solid curve in (a) and (b) is used for evaluating the signal yield. In (a) and (b), the vertical and 
horizontal lines denote ±5<r. Closed circles correspond to the data. The filled boxes show the MC 
signal distribution with arbitrary normalization. 

Assuming no correlation between them, all these uncertainties are combined in quadrature 
to give a total of 6.5%. 

Upper limits on the branching fractions at the 90% C.L. including these systematic un- 
certainties are calculated with the POLE program without conditioning The resulting; 
upper limits on the branching fractions are 

B(t~ -> e~K%) < 5.6 x lCT 8 
B{j~ -> yTKl) < 4.9 x 1(T 8 . 

SUMMARY 

In conclusion, we have searched for the lepton flavor violation decays r~ — > £~K§ [I = 
e or (j,) using data collected with the Belle detector at the KEKB e + e~ asymmetric-energy 
collider. We found no signal in either mode. The following upper limits on the branching 
fractions were obtained: B(t~ -> e~A|) < 5.6 x 10~ 8 and Bir~ —> \i~K^) < 4.9 x 1(T 8 
at the 90% confidence level and including systematic uncertainties. These results improve 
the search sensitivity by factors of 16 and 19, respectively, compared to the previous limits 
obtained by CLEO experiment. 
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